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1. Introduction 



Understanding nonperturbative aspects of Yang-Mills theory is one of the most important 
themes in modern theoretical physics. It has been playing the role of theoretical arena where 
many new ideas are tested. Also, in the attempts there have been obtained new analytical 
tools, each of which in turn always has opened a way to vast area of new applications. 

In the present paper, we investigate 3D Yang-Mills theory. This is certainly a simplified 
toy model for a realistic 4D theory, but still shares many of essential features with 4D YM 
theory, including infrared slavery, asymptotic freedom, color confinement and the existence 
of mass gap. Furthermore, the 3D theory has physically meaningful degrees of freedom, as 
opposed to more simplified 2D YM theory. 

Recently there have been remarkable developments in the understanding of 3D YM 
theory [1][2]. Among them is the study [1] based on the Hamiltonian formalism with the 
so-called Karabali-Nair (KN) variable [3] [4]. (For the analysis based on the Kogut-Susskind 
Hamiltonian see, e.g., [5].) 

The KN variable is a sort of Bars' corner variable [6], which is local and gauge- invariant. 
In their original papers [3], Karabali, Kim and Nair investigated the vacuum wave functional 
of 3D pure SU(N) YM theory, and determined its behaviors at the IR and UV regimes. 
They obtained the string tensions for finite iV's, which are in good agreement with lattice re- 
sults [8]. Recently, with a few Ansatze Leigh, Minic and Yelnikov [1] determined completely 
the Gaussian part of the vacuum wave functional in the 't Hooft limit, and obtained glue- 
ball mass spectrum. Their arguments are somehow intriguing, but the obtained spectrum 
exhibits excellent agreement with lattice data [8]. 

As will be reviewed briefly in Section 2, the most subtle part in the discussions of Leigh 
et al. is the renormalization of the Schrodinger equation, to remove divergences which always 
arise when two successive local variations act on the wave functional. This is usually carried 
out term by term on all possible terms in the wave functional, as was done in [1] . 

In the present paper, we propose a new algorithm to obtain the renormalized Schrodinger 
equation, which is certainly effective in the 't Hooft limit. We heavily use the fictitious 
holomorphic symmetry [3][1] existing in the framework based on the KN variable, and 
determine the possible local counterterms for the Lagrangian rewritten in terms of the KN 
variable. Then, in the 't Hooft limit the Schrodinger equation simply becomes the Hamilton- 
Jacobi equation for the resulting effective Lagrangian. 
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The equation thus obtained actually has a different form from the original one of Karabali 
and Nair, but for the Gaussian part of the ground state wave functional it gives the same 
result with that of Leigh et al. We believe that this gives a justification for their result. 
Since our equation is totally algebraic and easy to solve recursively in gauge potentials, one 
can easily determine higher-order corrections to the Gaussian part, as will be demonstrated 
in section 4. Inclusion of other fields and application to deconstructing 4D YM theory will 
be reported in our subsequent papers. 

The present paper is organized as follows. In section 2 we review the framework of 
Karabali, Kim and Nair [3], and stress the importance of local counterterms which may 
appear when the product of variations in the Hamiltonian are normal-ordered. In section 3 
we introduce a new algorithm to calculate such local counterterms, restarting from the bare 
Lagrangian rewritten in terms of the KN variable. We show that the form of local countert- 
erms are almost determined by a simple reasoning based on dimensions and holomorphic 
invariance. We expect that this is the simplest method to obtain them in the 't Hooft limit. 
We then write down the equation that determines wave functionals. In section 4 we show 
that our equation can be solved recursively in gauge potentials to give the ground state wave 
functional. We further show that the resulting Gaussian part coincides with that of Leigh 
et al. Higher-order corrections to the Gaussian part are also discussed. Section 5 is devoted 
to conclusions and outlook. A few Appendices are added where detailed explanations are 
provided for some formulas used in text. 

2. 3D Yang-Mills theory and Karabali-Nair variable 

We consider 3D Yang-Mills theory with gauge group G = SU(N). Spacetime coordinates 
are denoted by x = (x^) = (x,t) with x = (x l ,x 2 ), and gauge potentials by traceless 
antihermitian matrices belonging to su(N), 



Here the hermitian matrices t a (a = 1, ■ • ■ , N 2 — 1) are normalized as trt a t b = (1/2) 5 ab , 
so that f acd f bcd = N5 ab for [t a , t b ] = if abc t c . Note that for X = {l/i)t a X a and Y = 
(l/i)t a Y a , [X,Y] = (l/i)t a [X, Y] a with [X,Y] a = f^X b Y c , and [X,t b ] = t a (&dX) ab 



F^{x) = -t a F« v (x) = [d, + A^x), d v + A v {x)] = -F%{x) 



(2.1) 



(F^ = ^-^ + r bc 4<). 



(2.2) 
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with (adX) ab = f acb X c . We also denote X ■ Y = X a Y a = -2 tr(XY). We often write X ■ X 
as X 2 . 

We take the temporal gauge, A (x) = 0, after which there only exists the residual gauge 
symmetry, Ai(x,t) — > g~ l {x) At(x, t) g(x) (i = 1,2) with g(x) E G = SU(N). In order 
to define the Karabali-Nair (KN) variable from the remaining gauge potentials Ai(x,t), we 
first introduce complex coordinates z = x 1 + ix 2 , z = x 1 — ix 2 and often write them as 
x = (z, z). They give 

d = d z = ^{dx - id 2 ), d = d- z = h& + id 2 ), (2.3) 

and accordingly, we combine the gauge potentials as 

A(x) = A z (x) = - (Ai(z) - iA 2 {x)) , (2.4) 

A(x) = A- Z (x) = +iA 2 (x)) = -A\x). (2.5) 

There always exists a matrix-valued function V(x) G G c = SL(N, C) such that 

A = -dV ] A = V~ l dV (= -A*) . (2.6) 

One can easily show that V is unique up to the left multiplication of a holomorphic matrix: 

V(z t z t t)^h(z,t)V(z,z t t) {h(z,t) e SL(N,C)), (2-7) 

while the right multiplication of a time- independent SU(N) matrix corresponds to a residual 
gauge transformation under the temporal gauge: 

V(z, z, t) - V(z, z, t) g(z, z) {g(z, z) e SU(N) ) 

[U - g _1 Ag + g-'dg, A -> g _1 Ag + g~ l 5g) . (2.8) 

Theory thus must be invariant under 

z, t) -> h(z, t) V(z, z, t) g(z, z) (h(z, t) e SL(N, C), #(z, z) e SU(N)) . (2.9) 

The KN variable [3] is then defined by 

H{x) = V(x) V^(x). (2.10) 

This is a gauge-invariant, local variable, and transforms under the holomorphic transforma- 
tion (2.7) as 

H(z,z,t) -»• h{z,t)H(z,z,t)h(z,t) {h{z,t) = {h(z,t)) ] ). (2.11) 
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With the KN variable, the gauge potentials are written as 

A = V^dV + V~ 1 JV 1 

A = v~ l 3v, 

where the hermitian connection J[x) is defined by 1 

J = -dHH- 1 . 

This transforms under the holomorphic transformations (2.7) as 

J(z, z, t) -> h(z, t) J(z, z, t) h~ l (z, t) + h(z, t) dh~ l (z, t). 



(2.12) 
(2.13) 

(2.14) 

(2.15) 



Equation (2.6) gives the local change of variables from (A, A) to (V, V^), while eqs. (2.12)- 
(2.14) should be regarded as defining the local change of variables from (A, A) to (V,H), 
not to (V, J). This is because the latter viewpoint does not reflect the relation A = —A^. 

The Hamiltonian under the temporal gauge, Aq = 0, is then given formally by 

^ 2 5 5 2 



H 



Using the Gauss law operator 
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2 SA a (x) 6A a (x) g 



(ad£> 2 ) 



ab 



5A b (x) v ZJ 5A b (xY 
we have 6/6 A = (ad/^)- 1 (-adD z 5/5 A + J). Setting D = d + J, we have 

adD z = V" 1 adD V, adD s = V -1 add V, 

JA = ^Jj' TA = - V-(adS)-'adZ)A + (adZy-j, 

where V(x) = (V ab (x)) represents AdV(x), V(x) t b V~ 1 (x) = t a V ab (x). Thus, when acting 
on gauge- invariant wave functionals (X = 0), the Hamiltonian becomes 



(2.16) 
(2.17) 

(2.18) 
(2.19) 



// 



2 5J a 



(add)" 1 adD 06 -^r - %{dj a f 



5J b 



r 



(2.20) 



This bare Hamiltonian has only a formal sense because it always gives rise to a short- 
distance singularity (oc 5 2 (0)) when acting on a local functional of Jj(x), and thus needs to 
be normal-ordered. By writing the wave functional as 

2vriV 



*[J(x)] =exp 



L {g*N) 



(2.21) 



1 This J(x) differs from the original hermitian connection J(x) introduced by Karabali and Nair [3] by 
a multiplicative factor as J{x) = —(tt/N) J(x). 
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the Schrodinger equation becomes 

where m is the (dimensionful) 't Hooft coupling constant, m = g 2 N/2-K . The second term 
on the right-hand side is seemingly of order 1/N but cannot be neglected even in the 't 
Hooft limit (N — > oo with m fixed finite). This is because counterterms coming from planar 
diagrams are accompanied by the sum over loop index of order N so that the net result 
can be of order N° = (1/N) x N. Note that such additional terms should be local since 
they stem from short-distance singularities. Once such local counterterms are extracted, 2 
the second term can be neglected safely in the 't Hooft limit. 

One can guess the form of the local counterterms. We first note that finite local countert- 
erms in the Hamiltonian formalism will appear when two successive functional derivatives are 
applied to the same exponent (of the form Nti (something)) or when the functional deriva- 
tive acts on a regulator. Thus, the finite local counterterms in the Hamiltonian, i2i oc = 
-(g 2 /2) J [(5/ 5 A^ 2 ] loc , will take the form H loc ~ g 2 N J [Nf(A, d 3 A) + # (A, djA)6/6A?] , 3 
and will not include terms proportional to (5/ '6 A®) 2 . This can be thought of as a correc- 
tion to the original Lagrangian, whose form will be roughly estimated (using the relation 
Ei = (l/i)8/5Ai ~ {l/g 2 )Ai although this should be modified after inclusion of time- 
derivative terms) as 

Aoc ~ g 2 N 2 f(A, djA) + iNA$f?(A, djA) (2.23) 

with no terms proportional to A 2 . This expression will hold even when the theory is rewritten 
with the KN variable, so that we expect that the main part of the £i oc consists of dimension- 
two terms which does not include time derivatives (recall that g 2 N has mass dimension one) 
and dimension-three terms which is linear in the time derivative of a field. 

3. Effective Hamiltonian method 

Usually the renormalization in the Hamiltonian formalism is carried out with the following 

steps: 

2 Karabali and Nair claimed that there should arise an additional term in (2.22) as (1/2) J J ■ 8P/8J 

[3]. This agrees with our renormalizcd Hamiltonian (to be given later) in the IR limit. Further discussions 

on this point will be given at the end of section 3. 

3 The overall N comes from the loop-index sum, and the N in the parenthesis comes from the exponent. 
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1. Start from a given bare Lagrangian 

2. Obtain the classical Hamiltonian 



3. Solve the Schrodinger equation 

The third step needs renormalizing the wave functional term by term, which is actually a 
cumbersome procedure to carry out. 

Since we only need to extract finite local counterterms, the following algorithm to obtain 
the renormalized effective Hamiltonian turns out to be possible: 

1. Start from a given bare Lagrangian 

2. Find the effective Lagrangian including finite local counterterms 

3. Solve the Hamilton- Jacobi equation for the effective Hamiltonian 

As opposed to usual quantum field theories in four dimensions, the second step is actually 
manageable since 3D gauge theory is super renormalizable; local counterterms can be con- 
trolled in the 't Hooft limit simply by a symmetry argument and dimensional analysis, as 
we see below. 

We again start from the Lagrangian 



2g 2 

and consider the partition function 



(i»+L4 ,A] a ) 2 -(B a Y (3.1] 



l dA v\ JS YM [A 



Vol (gauge) 

We take the temporal gauge A (x) = 0, which leads to 



g^YMl^VJ . (3.2) 



Vol(res. gauge) 



Det^/H^oM^,] e i5YM[Aj 



Vol(res. gauge) 



with 



S Y m[A u A 2 ] = Jd 3 x^ (4 2 - B 2 ). 



(3.4) 
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With the change of variables (2.12)-(2.14), we have 

5 A = V~ l (5 J + [D, 5VV~ l ] ) V 

= V- 1 ( [D, -5HH- 1 ] + [D, 8VV- 1 ] ) V, (3.5) 
5A = V~ 1 d(6VV- 1 ) V, (3.6) 

where — 5V V~ x and — 5H H~ l are the Maurer-Cartan form, and we denote = —d^V V~ x . 
The metric over the functional space then becomes 

\\SASA\\ = - f tiSASA 

J X 

= fti (6VV 1 - (l/2)A- 1 BadD5HH- 1 ) ■ A [SVV~ l - (l/2)A- 1 dadD5HH~ 1 ) 

J X 

- (1/4) 6HH- 1 ■ adDadS A -1 adDadd 5HH' 1 ] , (3.7) 
where A is the Laplacian in the adjoint representation: 4 

A = - (adD add + add adD) . (3.8) 

The Haar measures of V and H are given by ||<5^|| 2 = — / tr(5VV -1 ) 2 and ||5iJ|| 2 = 
— J tr^SHH^ 1 ) 2 , respectively, and thus the functional measure becomes 5 

[dAdA] = [dVdH] ■ Det (-ad<9 ad£>) . (3.9) 

We also have 

A = V~ 1 (j- [D,u ])V, (3.10) 

X = V~ 1 (-duo) V. (3.11) 

Putting everything together, we obtain 

Z = [ \ dV( * H] - Det 1 / 2 ^ 2 ) • Det(-ad9ad J D) • e iSo ^ (3.12) 
J Vol (hoi. trsf) 



4 In the following ad 8 will often be abbreviated as 3. 

5 The measure [dVdH] (H = VV') is a natural extension of the area element dx A dy written in terms of 

,2 _ UI2. 



the complex coordinate z = x + iy = re and its absolute square r = \z\' 



J[dzd(r 2 )}---^ 1 - jTV) 



= / dx A dy 
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with 



S ° = 2^2 l d3x 



d A x 



-4(j-[D,uj }) -du G + A(dJ) 
-4u ■ Auo -Adu -j + 4 (dj) 



(3.13) 



One can easily see that So is invariant under the time-dependent holomorphic transformation 
(2.7)-(2.15), since tu transforms as u — > huj^h" 1 — hh" 1 , so that J — [-D,u; ] — > h(j — 

The Jacobian factor has the expected form e Arxtr ( somcthm g) ; anc [ thus the effective action 
with local counterterms will take the following form in the 't Hooft limit: 



S = S [uj , H] + Si oc [u , H] 



(3.14) 



wi 



th 6 



Si 



loc 



d 3 x £i oc (uj ,H, d„H) 



d x N 



a A 2 ) 



\3 r W 

'loc 



(3.15) 



Note that the finite local counterterms S'loc can depend also on u = —W~ x because it 
is invariant under time-independent (residual) gauge transformations. We have assigned 
dimensions to the fields as [uq] = 1, [H] = (thus [J] = 1), and demand S'loc to have the 
same symmetries with those of the original theory. That is, we require that Si oc be invariant 
under C, P (and thus T) as well as (time-dependent) holomorphic transformations. 

One can easily check that neither dimensionless nor dimension-one terms can contain 



nontrivial differential polynomials of the fields, so that we have C[® c = and C 



•(i) 

'loc 



0. 



As for dimension-two local counterterms, we find that £[H' oc dHH -1 ■ dHH" 1 satisfies 
the C and P invariance (see Appendix A). This term itself, however, is not invariant under 
holomorphic transformations (2.11) since each of the pieces dHH~ l and dHH' 1 does not 
transform homogeneously, 



dHH" 1 -> h dHH" 1 hT 1 + dh hT 1 , 
dHH- 1 -> h dHH' 1 h" 1 + Hdh h^H' 1 . 



(3.16) 
(3.17) 



6 There can actually be a dimension-four local counterterm, which corresponds to the countertcrm to 
(l/2<7 2 ) tv F 2 V . We assume that the renormalization of g 2 is already performed before the following discus- 
sions arc made. 
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An important point is that the inhomogeneous terms, dhh~ l and H dhh^H^ 1 , belong to 
the zero modes of the operator ad<9(= d) and ad-D, respectively, and thus, those pieces can 
be completed such as to transform homogeneously if we subtract proper zero modes from 
them: 

[dHH~\ ol = ld(dHH~ l ), [dHH-\ oX = -^ a dD(dHH^), (3.18) 

where the inverse of an operator is defined for its part without zero modes. We further note 
that there is a subtlety for such dimension-two terms (related to gauge mass terms) because 
inclusion of IR regulator can also give a finite (generically nonlocal) term. Expecting that 
the remover of unnecessary zero modes is supplied by a proper IR regulator, we set the 
dimension-two "local" counterterms to be 7 

£g> = [dHH-\ ol ■ [BHH-'U = JL^8(dHH- 1 ) .-^adDidHH- 1 ), (3.19) 

where (3 is a certain constant to be fixed later. By using the definition J = —dHH~ l and 
the formula 5 J = -ad£> (5H if" 1 ), c[H can be further rewritten as 

The dimension-three term is found to be 

= i — HH" 1 ■ dj [a is a certain constant to be fixed later). (3.21) 
2n 

This is actually invariant under C and P (and thus T) (see Appendix A), and is also 
invariant under (time-dependent) holomorphic transformations up to total derivatives: 

HH- l -dJ HH' 1 ■ dj -2dtr(hjh- 1 ) +2dtr(fih~ 1 H~ 1 dH). (3.22) 

Since HH" 1 can be written formally as HH~ l = — (ad-D) -1 ^/, the local counterterm could 
be expressed in terms of J as C{H = —i (a/2ir) (a,&D)~ 1 j ■ dj . 

We thus conclude that the "local" counterterms are given by 



'loc 



N 



g 2 NC[2 + C £ = /3 d" 1 5 J ■ (adD)" 1 8 J - ia (adD)" 1 J ■ dj. 

J 07T Z Z7T 

(3.23) 



7 Another way of saying this would be that we introduce an IR regulator which removes the zero modes 
of the Laplacian A = (l/2)(adDi9 + <9adZ)) from the Hilbert space. Note that in such Hilbert space, one 
can adopt the formal relations, d~ 1 d = 1 = 9<9 _1 and (adD) _1 ad£> = 1 = ad£> (ad!?) -1 . 
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Note that this has the expected form (2.23). At the present stage, this is simply a working 
hypothesis and should be justified only by directly calculating the functional determinants 
in (3.12), which we could not have done yet. 

The effective Lagrangian thus becomes 

C = Cq + C\oc 

2 / 2 N \ 

= — -wo ■ Au + — = Buj + ia — (adD^dj) ■ J 
9 ^9 27T ' 

^7T Z g 
The canonical momenta are calculated as 



o 2 N 2 - - 2 - 

+ P W B ^ Bj ■ (^D)~ l 9J + - {djf. (3.24) 



dC 4 2 - ■ 

vry = 7^— = -^Au + -dJ, (3.25) 
<9w # 2 g 2 

ttt=^ = — -Bu + ia — (adD^Bj, (3.26) 
g 2 2tt 

from which the (renormalized) effective Hamiltonian density 7i = i\y ■ ujq + ttj ■ J — C is 
obtained to be 

H = - 9 —d~ 1 (mj + a^- (adD^dj) ■ 1 

2 V Z7T J 

+ a 9 -^d- l dj- (in j + a^- (adD^Bj) - /3 9 -^ r B~ l 8 J ■ (adD)- 1 Bj - -Ad J)' 
Ait \ Ztt J 8tt z g A 

n 2 a 2 N / - - - - \ 2 - 

= -!Ldr 1 iTrj-X+aZ—- {{adD^dJ ■ d~ x T + B~ x Bj ■ tTijj - -{Bjf 



+ (a 2 -f3) 9 ^B~ 1 Bj-(adDy 1 Bj. (3.27) 



in 2 



Here 2 a (x) is the holomorphic Gauss-law operator that generates holomorphic transforma- 
tions: 

Z\x) = (adD)%ir»j)(x) = (adD)^-^-- (3.28) 

with (ad.D) afe = 5 ab d + j ach J c [x) . In Appendix B we show that the change of variables 
(2.12)-(2.14) leads to the following identity for a gauge-invariant functional P\J[x)\ which 
is written with original variables as P\J(x)\ = P[A(x) , A(x)]: 

X\x)P[J] = -d(y*(x) H^y 1 ), (3-29) 



11 



where V ab (x) is again a matrix element of AdV(x): V(x)t b V 1 (x) = t a V ab (x). 8 Thus, 



the Hamilton- Jacobi equation for the wave functional 

2ttN 



V[J(x)\ = exp 



(g 2 Ny 



P[J(x) 



(3.30) 



becomes 



2 



2 

m = 



Am? 5JJ 5 A 



a 



+ (a 2 - /?) — d- x dj ■ (&&D)- l dJ (m = g 2 N/27r) 



In Appendix B we further show that the operators 



Af 
AT 



J (adD^dJ - V-t = - jf ((adD)- 1 )^, y) Bj b {y) V a 



5 



SA c (x) 



5 



(3.31) 

(3.32) 
(3.33) 



count the number of d and D, respectively, in the gauge- invariant functional on which 
they act. That is, when P\J\ is expanded as a formal Laurent series in d and D, the 
operators Af and Af act linearly on each term in the expansion and return its orders in 
d and D, respectively. For example, on P[J] = J2 n c n J BJ ■ (a.dD) n Bj', they act as 
AfP[J] = E„ 2c n / dj ■ (adD) n 8j and AfP[J] = E„ (« + 2)c n / dj ■ (adD) n 8j. We have 
used the relation dj = [d, D\. Note that the operators and Af are parity-conjugates to 
each other and have a definite meaning that is invariant under holomorphic transformations. 

One can see that the last term in (3.31) leads to an IR divergence in the gauge invariant 
sector since dj represents the field-strength, dj = (—i/2) VF^V' 1 . We thus set (3 = a 2 
and drop the term in the following discussions. The final form of the Hamiltonian density 
and the corresponding H-J equation thus become 9 

,2 „2 AT , 2 



H = - y <T 1 in j ■ 1 + a ^ ([&dD)~ l dJ ■ d- l T + d~ x dj ■ m a 



g 1 



(dj) 2 , (3.34) 



8 Notc that V = (V ab ) is an orthogonal matrix, (VX) ■ (VY) = X ■ Y. 

9 Another possibility for dimension-three local counterterm other than (3.21) could be to take the WZ-likc 
term : 

C^l = const tr (HH- 1 [dH if" 1 , dH if" 1 ]) . 

This is actually invariant under (time-dependent) holomorphic transformations because infinitesimal holo- 
morphic transformations (actually any local changes of variable) give total derivatives. One can easily check 
that this is also invariant under C and P (and thus T) (see Appendix A). We expect that it yields the same 
Hamiltonian with (3.34) when rewritten in terms of J together with a proper IR cutoff corresponding to 
this local counterterm. 
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l E ^Sf P+ - iMP .J[A ? ^. V f A + (Bjf] ( m = ^/ar). (3.35) 

We conclude this section with a comment on the value of a. The preceding discussions do 
not fix the actual value of a. For example, the removal of the IR divergence can be performed 
for any value of a (when (3 = a 2 ). Although a should be determined by directly evaluating 
the functional determinants in (3.12), this can be fixed by comparing our local counterterms 
with the finite corrections obtained by Karabali and Nair in their original Hamiltonian 
formalism [3]. To see this, we first notice that in spin-0 wave functionals, 5 always appears 
with the combination dj = [d, D] or in the form of A = (l/2)(adL><9 + <9ad£>). Thus, 
our operator Af can also be interpreted as counting the number of D (i.e. Af = A/"), or the 
number of J if we neglect d in D = d + J in the IR limit; Af = Af ~ J J - 5/5 J (IR). 
Then (3.35) in the IR regime agrees with the one originally obtained by Karabali and Nair 
using the heat kernel regularization [3] if we take a — 1. We set a = (3 = 1 in the following 
argument. 



4. Ground state wave functional 



In this section, we solve the H-J equation for the ground state wave functional ^[^(cc)] 



exp 



{2nN/(g 2 Ny)P[J(x)}} : 



= -MP + -NP 

4 4 



1 5P „5P 

^5j- v M + ^ 



(4.1) 



4.1. Preparations 

We expand the functional P[J] as 

oo 



1 



n=2 



(4.2) 



Here P n [J] (n > 2) is a holomorphic invariant functional whose order in J is at least n. 
For example, 

HJ] = jaj.K{^)3j = j 

is the Gaussian part augmented with higher order terms in J such that it becomes holo- 
morphic invariant. This can be rewritten in terms of the original gauge potentials as 



(^|) BJ + 0(J 3 ] 



(4.3) 



P 2 [J] = J [D- z ,D z ]-K^yD- z ,D z ] = P 2 [A,A], 



(4.4) 
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where A is the Laplacian with respect to the original gauge potentials: 



A = - (acID z ad£> z - + adD z ad£> z ) = (AdV)- 1 A AdV. 



(4.5) 



Since P n \J\ is holomorphic invariant, it must be constructed only from dj and A = 
(1/2) (<9D + -D<9), and thus, when 5j5J acts on P n [j7"], we only need to consider two cases: 
one is when it hits on dj in P n \J\ and the other is on A. Both cases can be easily dealt 
with by using the following formulas: 

5 



5J 



dJ-X 



-dx 



-a s 



5dj 



dJ-X 



j ^■y = \([d^Y]-[XM]), 



(4.6) 
(4.7) 



the latter of which can be derived as 



5j J AX -Y = 1 J(ad6J -d + d-ad5J)X -Y = J 5J ■ ^([dX,Y] - [X,dY]). (4.8) 

Note that the former case (when 5j5J hits on dj) necessarily reduces the order of the 
functional in J by one, while the latter (when 5j5J hits on A) does not change the order. 
This observation leads to the formula: 

~5P n 



SPn 


\5P n 




\SP n i 


5J 


V5J. 


+ 

n-1 


16 j\ 



-o 



+ 



5dJ L 5 J \ n 



(4.9) 



Similar analysis can be made for V5/5A acting on P\3\ = P[A,A]. We then obtain the 
following formulas: 

v jxj Bj ' x = v JlS ][ - D ~ zlD ^' jc = MX (~ adD Wj J fBj ' x )> (410) 

V-LJ AX-Y = V-LJ AX-Y = ^([ a dDX,Y] - [X,adDY]), (4.11) 

which are derived as 

5 A J [D Z7 D Z ] -X = j [5A,D Z ] - X = j 5A- [D z ,-X], (4.12) 

5 A J AX ■ Y = ^ J ( &d5A &dD z + &dD z &d5A )X ■ Y 

= J 8A-^([adD z X,Y] - [X,adD,y]). (4.13) 

We also have a formula similar to (4.9), 
.5P n 



V- 



5A 



v 5 ^ 

5 A J n-l 



+ 



5A\n 



SdJ 



5 A Jn 



(4.14) 
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For example, on P 2 [J] = J2( K n/m 2n )O n with O n = j dj ■ A n dJ, the operators 8/8 J and 
V 5/5 A act as (setting L = A/m 2 ): 



5P2 

5 J ' 
with 



8_Pz 
.5 J 

8JV\ 
5J 



+ 



5Po 



5J. 
-2dK{L) 8 J, 



5JV\ 
5J 



oo n— 1 



K 



EEd^ Am ^ An ~ m ^L ( 4 - 15 ) 



2 — ' — ' m 

n=0 m=0 



and 



.5Po 



5A 



with 



V 



. 5A 
5_h 
5A 



V 



.5P2 

5A 



2adDK(L) dj, 



.6P2 
5A 



00 n—l 



EE^[ adDAm ^ An ~ m "^]- 

(4.16) 



2 L — ' *■ — ' m 

n=0 m=0 



4.2. Gaussian part 



After making rather lengthy preparations presented above, we are now in a position to 
calculate P n \J\ {n > 2). First, by using (4.15) and (4.16), the holomorphic-invariant 0(J 2 ) 
part in (4.1) is expressed as 



4 4 4m 2 



5Po 



5J 



V 



6Ah 



- / {djf + 0(J 3 



= -^MP 2+ l -MP 2 - J 8 J ■ LK 2 (L) 8j- J (8 J) 2 + 0(J 3 ). 

For K(L) = J2n K nL n with L = A/m 2 , Af acts as 10 

N J 8j -K{L)8j =AfJ2 K " J Bj - (^) n ^ = E( n + 2 )^/ a ^ 



(4.17) 



— ) M 



&7- [L— + 2)K(L)dJ. 



The action of A/" is the same with that of M for functionals of spin zero. 

Af j 8 J ■ K{L) 8j= j 8j- (L-£- + 2)K(L) 8 J. 



(4.18) 



(4.19) 



10 This gives a justification for the rule given by Leigh et al. (see Appendix A of [1]), where (4.19) is 
assumed to hold as a consequence of complicated normal ordcrings in subject to successive variations in the 
Schrodingcr equation. Note that their operator J J-8/5J is replaced by (N +AT) /2 in our framework, where 
we assume that any necessary normal orderings are already taken care of by adding local counterterms to 
Lagrangian. 
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Thus, collecting the coefficients of J 2 terms in (4.17), we obtain the differential equation 
for the kernel K(L), 

d 



(^tijiM-^D-m 



(4.20) 



which certainly reproduces the result of Leigh et al. [1]. 

Once the equation (4.20) is obtained, one can derive all the results in [1] with the same 
logic. First, the equation has a normalizable solution 

1 J 2 (4v / I) 



K(L) 



VL J!(4x/I) : 



(4.21) 



which has the following asymptotic form: 

K(L) ■ 1 + (2/3)L 



(L -> 0: IR) 
(L -> -oo: UV) . 



Thus, the ground-state wave functional behaves in the IR and UV regions as 



(4.22) 



*ir = exp 
= exp 

^uv = exp 
= exp 



2ttN 

7T 



2g*N 

2 -Joj 



(Bj(x)f + - 



(4.23) 



-an 



Bj 



l 



F 12 (x) 



-an 



F 12 (x) 



(4.24) 



The UV wave functional reproduces that of massless free gluons. The vev of a large Wilson 
loop, on the other hand, can be evaluated with the use of ^ir as 

(trPe"^ ~ J[dA{x)dA{x)\ trPe"^- |^ IR) | 2 

= J [dA(x)dA(x)] tr Pe~ §A ■ e~ {n/9iN) U F *&»* + •". (4.25) 

This is nothing but the Wilson loop average of 2D YM with 2D gauge coupling g\ = 
g A N/2ir = mg 2 . As was first pointed out by Karabali et al. [3], the string tension of 3D 
pure YM theory is then calculated from the known value of the string tension of 2D YM 
theory [7] as a = g\ (N 2 - I) /AN = g 4 (N 2 - l)/8vr, i.e. 



g 2 N 



57T 



— m (N — > oo). 



(4.26) 



The computation of 3D glueball masses could be carried out in the same way as in [1]. 



16 



4.3. Higher-order corrections 

As for the 0(JT 3 ) part, we have the equation 



= l -NPz + l -NPz 



1 

2 



~MP 3 + l -UP, 





\SP 2 - 




/( 


lsj\ 


1 



5A m 2 5A 



V 



5P2 

5A 



5Po 1 5P« 



i 15J m 2 8J 



oo n—1 



(-BK(L) 5 J) ■(££-£ [adD A"W, A— ^ + ^adD 



n=0 m=0 

OO 77 — 1 



58 J J 



71=0 777 =0 

OO OO / — 1 

ZEE 



^ ^ fc=0 1=0 771=0 



where we have introduced the operators 



K k K t , 



Cfe,m,i-m-l ) + 0(^7" 

(4.27) 



(9 



k,Lm 



J A k dJ ■ (8 [adDA l 8J, A m 8j] + adD [A m dJ, dA l dj] ) (fc, Z, m > 0). (4.28) 
They are totally symmetric with respect to their indices up to 0(J~ 4 ). Actually, 
O kjlim = J A k Bj ■ (8[&dDA l 8j, A m Bj] + adD [A m 8j, dA l Bj] ) 

= J A k 8j ■ ([&dDA l 8j,8A m 8j] + [&dDA m 8j,8A l 8j]) +0(J 4 ) 
A k 8j ■ ([adDA m Bj,BA l Bj] + [adDA l dJ, 8A m 8j] ) + 0(J A ) 

A k Bj ■ (3[adDA m dJ, A l Bj] + adD [A l Bj, 8A m dj] ) + 0(J 4 ) 
= O k , m>l + 0(J A ) : (4.29) 
O kjlim = J A k 8j ■ (8 [DA l 8j, A m 8j] + adD [A m 8j, 8 A 1 8 J] ) 

= J (adDA l Bj ■ [BA k Bj, A m 8j] + 8 A 1 8 J ■ [-&dDA k 8j, A m Bj] s j 
A 1 8 J ■ (adD [A m 8j, 8A k 8j] + 8[adDA k 8j, A m 8j] ) 
CW (4.30) 
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Equations (4.27) and (4.28) motivate us to set P 3 to be in the following form : 

Ps = m 2(fci+m+l) M,m (4-31) 

with totally symmetric coefficients d k j )Tn . Then, up to 0(j7" 4 ) terms, each part in (4.27) is 
rewritten into the following form: 

(k + I + m + 4) 



klm 

K(L) 9J ' L 6dJ = 2^2^ m 2(n +k+ l +m+ 2) I A ^ ' 



S m 2(n+fc+i+m+2) (°*+n+l,I,m + Cfc.Z+n+l.m + Cfc,/, m+ n+l) 

^ m 2(fc+/+m+l) + dk,l-l-n,m + dk,l, m -l-n)^k,l,m, 

n k.l.m 



(4.33) 



\ ^ \ ^ \ ^ K r K n _ \ ^ KkKi +m+ i 

m 2(r+n) U r,n-l-m,m - m 2(fc+;+ m+ l) U M>m 

r n m=0 klm 



1 ^ (K k K l+m+1 + KiK m+k+1 + K m K k+l+ i) 

k,l,m 



o 



k,l,m- 



3 m 2(fc+«+m+l) 

(4.34) 

Collecting the coefficients of O k ,i, m , we obtain the recursion equations for c4,z,m : 
= (k + I + m + 4)cZ M>rn - 2 (Y^ K nd k —l—n, I, m + Y, K ndk,l — l—n, m , m—l—n I 

n n n 

— 7^(KkKi +m+ i + KiK m+k+ i + K m K k+ i +1 ), (4.35) 

which can actually be solved recursively for the increasing ordering in the sum of indices, 
k + I + m, with the initial condition rf ,o,o = (l/4)-K"i = 1/6. One can easily check that this 
recursion equations are equivalent to the following differential equation for the generating 
function D(Li, L 2 , L 3 ) = T,k,i,m d k,i,mLi k U L 3 m : 

r\ r\ r\ 

L x (— - 2K{L x Yj + L 2 (— - 2#(L 2 )) + L 3 (— - 2K(L 3 )) + 4] D(L U L 2 , L 3 ) 
= 3l A(Ll) L 2 -L 3 +K{L2) L 3 -L X +A(La) L,-L 2 )■ (436) 
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Calculation of 0{J A ) part can also be performed in a similar way. 

We conclude this section with a comment on the glueball mass calculations carried out by 
Leigh et al. [1] , who set the Ansatz that the vacuum wave functional is Gaussian. Although 
the obtained results are in excellent agreement with lattice data [8], we have seen above 
that there exist nonvanishing corrections to the Gaussian part, so that there must be some 
explanation why the Gaussian Ansatz gives a good approximation. We still do not have a 
satisfying answer for that yet, but we expect that higher-order corrections will take the form 
of integrating a rational function of Bessel functions, which will get suppressed strongly due 
to the cancellation in oscillating integrals. 

5. Conclusions and outlook 

In this paper we proposed a new method to analyze 3D Yang-Mills theory in the 't Hooft 
limit based on the effective Lagrangian and holomorphic symmetry. We wrote down the 
equation which determines the wave functional, and showed that it can be solved recursively 
in J . We gave sample calculations for the Gaussian and 0(J^ 3 ) parts of the ground state 
wave functional, and showed that the Gaussian part reproduces the result obtained by Leigh 
et al. [1]. 

There are many things that should be carried out. First of all, we need to elaborate 
the analysis itself, especially to further clarify the origin of the local counterterms and to 
determine the coefficient a in (3.35) by directly evaluating the functional determinants in 
(3.12). 

Other possible directions are to consider finite size effects in the present framework, and 
also to include many other fields including scalars and fermions in various representations. 
Among them, the inclusion of adjoint scalars should be particularly interesting, since it can 
be used to construct 4D YM theory via deconstruction [9]. Investigations along the above 
line is now in progress and will be reported in our future communication [10]. 
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A. J pc 

We summarize the quantum numbers assigned to fields. 
J 

A field <&(z, z) is said to have spin J when it transforms as 

&(z,z) =e ije §(e id Zj e- ie z). (A.l) 
This definition leads to the following assignment: 





A A d d H J 


dj D A 


J (spin) 


+ 1 -1 +1 -1 +1 


0+10 



E: 

Parity transformation is defined as (x 1 ,^ 2 ) (x 1 , —x 2 ) or 

x—(z,z,t) — > xp = (z,z,t), (A. 2) 

under which the original gauge potentials A(x) and A{x) transform as 

A{x) Z> A{x P ) (A.3) 

A{x) Z> A{x P ). (A.4) 

This gives 

V(x) ^ V^-\x P ) H{x) 4 H~ x {x P ) (A.5) 

and 

cu (x) ^ + (H- x {u Q + HH- x )H)(x P ) (A.6) 

J(x) ^ + {H- l dH){xp) (A.7) 

dj{x) £ - [H~ x dJ H)(x P ), (A.8) 
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which show that the matrix HH l dJ transforms as 

HR- x dJ{x) ^ + (H^HH^dJ H)(x P ). (A.9) 

C 

Charge conjugation is defined for the original gauge potentials A(x) and A(x) as 

A,(x) Z - A J(x) (A.10) 
A,(x) Z - A T Jx). (A.ll) 



This gives 



and 



V(x) {V T )-\x) => H{x) $ {H T y\x) (A.12) 



u (x) Z - u^(x) (A.13) 

J(x) Z - J T (x) (A.14) 

Bj(x) Z - Bj T (x), (A.15) 

which show that the matrix HH~ x Bj transforms as 

HH^Bj(x) Z (Bj HH~ 1 ) T (x). (A.16) 



B. Holomorphic Gauss- law operator X a (x 

We collect properties of the holomorphic Gauss-law operator defined in (3.28), 

where (ad.D) ab = 5 ab d + f acb J c (x). The following statements actually hold for arbitrary 
finite N. 

We first show that under the change of variables (2.12)-(2.14) the variations transform 



as 11 



l 1 = V~ 1 ( a dd)- 1 (-l+ pi, (B.3) 



SA V Su 



11 In text, ad d is usually abbreviated as d. In the following, the variation with respect to H is always 
rewritten in terms of J by using the relation 6 J = —a,dD (SH H~ l ). 
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where 5cu(x) = — 5V(x) V 1 (x), and V(x) = (V a6 (a?)) is the matrix representing AdV(x); 
V(x)t b V- 1 (x) =t a V a \x). 

[proof] 

In general, when a set of fields {$/(cc)} are (locally) transformed to another set as 
6$i(x) = V u 5<f>j(x), V = {V u ) = dn-i^x) d h ■ ■ ■ d in , (B.4) 

their variations are transformed as 12 

(^)utAt^ (B.5) 



/(pt-i )/J (^^_J_ ) (b.6) 

J V 



5(f>i{x) v ' J 5<f>j(x 
^TT^= l^- l ) LJ {x^ 

where pt = ^ n (-l) n c\ • • • d in d£... in (x). We apply this formula to the sets {$/} = {A a , A a } 
and {0 7 } = {J a ,ou a } with <5cj a = -(iyr 1 )" Equations (3.5) and (3.6) give 

13 = Pv) = ! * > (B.7) 

\ V- 1 / \ -add / 

so that we obtain 

*=( 1 °)f V °), (B.8) 
\ adD add / \ V J 

- ( V " ° V - 1 ° V (B.9, 

\ V- 1 / \ -(add^adZ) (add)" 1 J 

Here we have used that V is orthogonal, (V T ) _1 = V, and ad-D and add are anti-self-adjoint, 
(ad.D)t = — a.dD and (add)''" = —add. Substituting these into (B.6) and noting that 

&dD z = V" 1 ad£> V, ad£> 2 = V" 1 5 V, (B.10) 

we obtain (B.2) and (B.3). [Q.E.D.] 

Note that the original Gauss-law operator that generates gauge transformations, 
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This can be proved by noting that the following holds for an arbitrary functional G[$]: 
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can be rewritten with new variables as 



i\x) = (V-^r'^y {5u = -5VV-1). (B.12) 

Thus, a functional is gauge-invariant when and only when it does not depend on V. Note 
also that (B.3) implies that 

1 = -add V-r + — = -add V-r + VI. (B.13) 

5A 5uj 5A v ' 



TP = -B(V 5 -^]. (B.14) 



Thus, for a gauge-invariant functional P\3\ = P[A, A] (with IP = 0), we have 

5P\ 
Si) 

This is the formula (3.29) used in text. 

We now investigate properties of the operator defined in (3.32), 

N=-J (^Dy'dJ-V-t = - jf {( a dD)- 1 ) ab (x,y)dJ b (y) ^J^y (B.15) 

when it acts on a holomorphic invariant functional P\J\ = P[A,A]. In general, P[J] can 
be expressed as a formal Laurent series of d and D. Since J\ is a first-order functional 
derivative, it obeys the chain rule as well as the linearity, and thus we only need to know 
the behavior of Af when acting on a polynomial of d and D of the same order, like 

P[J] = J dJ-(zdD) n dJ {n E Z). (B.16) 

In this example, the orders in 5 and D are 2 and n+2, respectively (recall that dj = [d, D]). 

We prove that M counts the order in d of each term in the Laurent series. We first note 
that when rewriting P[J~] as a functional of the original fields A and A (i.e. P\J\ = P[A, A]), 
d in P\J\ is replaced by ad-D 2 = ad(<9 + A) in P[A, A]. This is the only possible way for A 
to appear in a gauge-invariant way. For example, the polynomial (B.16) is rewritten as 

P[A, A} = J [D g , D z ] ■ (adD z ) n [D- z , D z }. (B.17) 

Due to the fact that the operator M is a first-order functional derivative with respect to A 
and also that it obeys the chain rule, the behavior of the operator M acting on a functional 
of P[J] = P[A, A] can be known once we understand how the operator acts on a functional 
of simpler form 

P'[A,A] = J X-[B + A,Y] 

( = P'[J] = J X ■ BY with X = V~ Y X and Y = V^Y^j (B.18) 
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with X and Y being independent of A. We emphasize again that this simplified functional 
is a representative of the part on which the variation 5/5A acts. It is thus sufficient to show 
that MP' = +P', which will be proved as follows. 

We first notice that 

MP' = - J (adD^Bj ■ V-L J X- [B + A,Y] 

= - J v- 1 {&&Dy l dj ■ -L J X ■ [B + A,Y] 

= + J (adD,)- 1 [D„D S ] ■ [Y,X]. (B.19) 

One could evaluate this by formally deforming the last expression as J (ad-D^) -1 adD z (D s ) ■ 
[y, X] = f D ? ■ [y, X] — J [D 2 , y] ■ X = P'. In order to justify such formal manipulations, 
especially to give grounds to the second expression, we make a few preparations in the 
following. First we extend the Lie algebra su(iV) to u(iV) and introduce the set of u(N)- 
valued differential operators: 

iS = {u(A r )-valued differential operator}. (B.20) 

An element T> e S generically has the following form: 

V = -t a V a (a = 0,1,--- ,N 2 - 1; t° = ^=l N ), (B.21) 

where T> a is a differential operator. For example, the covariant derivatives D z and D 2 belong 
to S and have the form 

D z = l N d + -t a A a , D- z = l N d + -t a A a (a = l,--- ,N 2 -1), (B.22) 

i i 

where d and A a are differential operators with components d = ((cc Idly) = d5 2 (x — t/)) 
and A a = ((x\A a \y^ = A a (x)S 2 (x — t/)), respectively. We introduce the trace for S as 

™ s {V€S) ' (R23) 

which satisfies the relation TtViV 2 = TrP 2 ^i for T>%, T> 2 € S. The action adD for D G S 
is then regarded as belonging to End(iS), and similarly (adP) -1 G End(5) when its inverse 
exists in any formal sense (which we always assume in the presence of a proper IR regulator). 

Now (B.19) can be interpreted as taking a trace in S and can be rewritten in the following 
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way: 



MP' =-2 J tr ((acLD,)- 1 [D„D S ]) [Y,X] = -2 Tr ((adD,)" 1 &dD z [y, X] 

= -2 Tr [y,l] = -2 Tr X [D 2 ,y] = ^ X ■ [D S ,Y] = j X ■ dY 
= P\ (B.24) 

which completes the proof of our claim. 

Similarly, we can prove that M = J d~ x dj ■ 5/6 J counts the order in D of each term in 
the Laurent series. It is again sufficient to show that MP" = +P" for a functional 

P"[J] = j X ■ [D,Y] = ! X ■ [d + J,Y] ((X,y): ^-independent) . (B.25) 

J X J X 

This can be performed as follows: 

MP" = j B ~ lBj 'Jj j x x ' [ d + J> Y ] 

= j d~ l [d,D] ■ [Y,X] 

= -2 Tr ((adS)" 1 add (£))) [Y,X] 
= -2 Tr D [y X] = -2 Tr X [D, Y] 

= P" . (B.26) 
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